Dirac Fermions in Graphene Nanodisk and Graphene Corner: 
Texture of Vortices with Unusual Winding Number 
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We analyze the zero-energy sector of the trigonal zigzag nanodisk and corner based on the Dirac 
theory of graphene. The zero-energy states are shown to be indexed by the edge momentum and 
grouped according to the irreducible representation of the trigonal symmetry group Csv Wave 
functions are explicitly constructed as holomorphic or antiholomorphic functions around the K or 
K' point. Each zero-energy mode is a chiral edge mode. We find a texture of magnetic vortices. 
It is intriguing that a vortex with the winding number 2 emerges in the state belonging to the E 
(doublet) representation. The realization of such a vortex is vary rare. 



Graphene nanostructure[l| has opened a field of 
carbon-based nanoelectronics and spintronics alterna- 
tive of silicon or GaAs. Carbon is a common mate- 
rial and ecological. Large spin relaxation length is ideal 
for spintronics [3. The basic graphene derivatives are 
nanoribbons[3|-|5| and nanodisks[6 10). They correspond 
to quantum wires and quantum dots, respectively. Re- 
garding a nanodisk as a quantum dot with internal de- 
grees of freedom, similar but richer physics and applica- 
tions are anticipated[I3|- There are many type of nan- 
odisks, among which the trigonal zigzag nanodisk [FigfT] 
is prominent in its electronic and magnetic properties 
owing to the zero-energy sector The main feature 
is that it becomes a quasiferromagnet in the presence 
of Coulomb interactions [6]. Recently trigonal nanodisks 
have been experimentally obtained by way of the Ni 
etching of a graphene sheet An experimental real- 
ization will undoubtedly accelerate further experimental 
and theoretical studies on graphene nanodisks. 

In graphene, the physics of electrons near the Fermi en- 
ergy is described by the massless two-component Dirac 
equation or the Weyl equation [pi - HoT ] . Nanoribbons have 
been successfully analyzed based on the Weyl equation 
but this is not yet the case with respect to nanodisks. 
The Dirac theory of graphene nanodisks must be indis- 
pensable to explore deeper physics and promote further 
researches. 

The carbon atoms form a honeycomb lattice in 
graphene. We take the basis vectors ai = (1,0) a and 
&2 = (l/2,\/3/2)a with a the lattice constant (a « 
2.46A). The honeycomb lattice is bipartite and has two 
different atoms per primitive cell, which we call the A 
and B sites [FigfT]. The Brillouin zone is a hexagon in 
the reciprocal lattice with opposite sides identified. We 
start with the tight-binding model (TBM) only with the 
nearest neighbor hopping t. We ignore the spin degree of 
freedom in most parts in what follows. 

The band structure is such that the Fermi point is 
reached by six corners of the first Brillouin zone, among 
which there are only two inequivalent points. We call 
them the K and K' points. The dispersion relation is 
linear around them, E K (k) ~ fi/i>F|k — K T | with r = ±, 



FIG. 1: (Color on line) Trigonal zigzag graphene nanodisk. 
(a) The nanodisk size is defined by N — N\, cn — 1 with iVben the 
number of benzenes on one side of the trigon. Here, N = 5. 
The A sites are indicated by red dots. The electron density 
is found to be localized along the edges, (b) A vortex texture 
emerges in the real-space Berry connection. In this example, 
the winding number is 2 for the vortex at the center of mass 
while it is 1 for all others. 



where vf = V3ta/2h is the Fermi velocity and K T = 
a" 1 (t27t/3,27t/V3) for the K point (K+) and the K' 
point (K-). 

The dispersion relation near the K and K' points is 
that of 'relativistic' Dirac fermions. Indeed, the TBM 
yields the quantum-mechanical Hamiltonian[13l4l5j. 





7~kx ~ I - ^kfj 







(1) 



where we have introduced the reduced wave number by 
k = \i—K T . The Hamiltonian acts on the two-component 
envelope function, $ r = (4>\,4>r). Each Hamiltonian 
describes the two-component massless Dirac fermion or 
the Weyl fermion. The Weyl equations read 



ihdt$ T { x ) — vya ■ p T <I> r (x), 



(2) 



where p T = h(Tk Xl k y ) = —ih(Td x ,d y ). The wave func- 
tion is given by i/'s( x ) = e m 4>s{x)- 

The symmetries are as follows. We note that Hk> = 
cfyH^cTy and a z H T a z = —H T , where a y and a z are the 
generators of the mirror symmetry and the electron-hole 
symmetry, respectively. 

In terms of the complex variable, the Weyl equation 
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FIG. 2: (Color on line) (a) Trigonal corner of graphene. The 
electron density is found to be localized along the edges, (b) 
Real-space Berry connection for the trigonal corner. A series 
of vortices are found to be present. 

reads 

d z *cf>l(x) = ie<f%{x), d z <j>l{x) = U</%(x), (3a) 
d,<f%(x) = -fe^f (*), d z *</>%'(x) = -ie<f%(x), (3b) 

with e ~ E/2Hvy- The envelope functions are holo- 
morphic or antiholomorphic for the zero-energy state 
(E = 0). 

Dirac electrons on zigzag edge: We analyze a graphene 
sheet placed in the upper half plane (y > 0) with the 
edge at y = 0. Translational invariance in the x direction 
dictates the envelope function is of the form 4>g{x,y) — 
e 1 ** 1 fg (y). Due to the analyticity requirement wc obtain 

(t>f( x ) = Cfe lk{x+ty \ <j%{x) = c£e fk{x ~ lv \ (4a) 
cj)f(x) = Cfe fk{x - iy \ tf$\x) = Cf^ x+iv \ (4b) 

with Cg being normalization constants. Hereafter we 
ignore such normalization constants. 

According to the TBM result, there are no electrons 
in the B site on edges. Hence we require 0g (y = 0) = 0. 
By avoiding divergence at y — > oo, the resultant envelope 
functions are found to be (f)\{x) — e lkz for k > and 
ift^(x) — e lkz for k < 0, with all other components 
being zero. 

The wave number is a continuous parameter for an 
infinitely long graphene edge. According to the TBM 
result, the flat band emerges for 

— 7r < ak < —2ir/3 and 27r/3 < ak < tt, (5) 

around the K' and K points, respectively. The boundary 
points ak — — 7r and ak — n are to be identified since 
they represent the same point in the Brillouin zone. 

Dirac electrons in trigonal corner: We apply the above 
analysis to the study of the envelope functions for elec- 
trons in the zero-energy sector of the zigzag trigonal cor- 
ner [Fig^a)], which is an infinite region surrounded by 
the x axis and the line with the angle argz = 7r/3. They 
are holomorphic (antiholomorphic) around the K (K') 



point. Here we discuss envelope functions around the K 
point. We start with the solution <\>\ (z) = Ce %kz for the 
upper half plane. We rotate this by the angle 7r/3, which 
presents us with the solution 2 (z) = £7 e lfc2:cx p[2W3] f or 
another half plane. The trigonal corner is given by the 
overlap region of these two half planes, which is described 
by a linear combination of these two functions with an 
appropriate coefficient. It is to be fixed by imposing the 
boundary condition: Since the top of the corner is located 
at z = 0, where there is no atom, we impose </>a(0) = 0. 
The resultant function is (f)\ (x) = </> (z) with 

(z) = e lkz - e lkz cx p[ 27 ™/3] _ ( 6 ) 

The phase shift is 7r at the corner. The envelope function 
around the K' point (k < 0) is given by </>a'( x ) = ^ ( z *)- 
We calculate the real-space Berry connection, Ai (x) = 
—i{(j)\di\4>). It exhibits a series of vortices [FigJ^b)], 
where the wave function vanishes. It has an infinite 
number of zero points at z n — (27m/v / 3fc)e 7rl / 6 , n = 
1, 2, 3, • • • , around which it is expanded as (f> (z) — k(z — 

Zn)- 

Dirac electrons in trigonal nanodisk: Our main pur- 
pose is to apply the above result to the analysis of the 
zero-energy sector of the trigonal zigzag nanodisk [Figfl] . 
The envelope function of the trigonal zigzag nanodisk can 
be constructed by making a linear combination of enve- 
lope functions for three trigonal corners. We consider the 
trigonal region whose corners are located at z\ — (L,0), 
z 2 = (— L,0), z 3 = (0, "\/3£)- As the boundary condi- 
tions we impose 4> ( z i) = & i z 2) — <t> (23) = 0. The enve- 
lope function is obtained around the K point (k > 0) as 

0K( X ) = 0( z ) w i t h 

<f)(z) = e ikz — e ikL e ik{z-L) cxp[-27ri/3] 
ikL^ik(z+L) cxp[27ri/3] 

The envelope function around the K' point (k < 0) is 
given by 0a' (x) = (f>(z*). Note that 0g(x) = identi- 
cally, as is consistent with the TBM result [6]. 

The wave number is quantized for a finite edge such as 
in the trigonal nanodisk. We focus on the wave function 
-0a ( x ) a t one of the A sites on an edge. For definiteness 
let us take it on the 2- axis. We investigate the phase 
shift between two points (x, 0) and {x + ma, 0), 

27T771 

6 T (x, m) = r— h arg0 (x + ma) — arg</> (x) , (8) 

o 

with ([7]). There are N links along one edge of the size- 
N nanodisk, for which we obtain precisely r (a/2, N) = 
Nak. On the other hand, the phase shift is tt at the 
corner. The total phase shift is 3Nak + 3ir, when we 
encircle the nanodisk once. This phase shift agrees with 
the TBM result. By requiring the single- valueness of the 
wave function, and taking into account the allowed region 
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FIG. 3: Vortex textures in the real-space Berry connection for 
the state \k") in the nanodisk with N = 7. The representation 
is indicated in the parenthesis. There are n vortices along the 
y-axis in \k"). A vortex appears at the center of mass for the 
15- mode It is interesting that the winding number is 2 

in the state \k„)- 

of the wave number ([S]), we find that the wave number is 
quantized as 

ak n = ± [(2n+ l)/3N + 2/3]n, < n < (N - l)/2. 

(9) 

When N is even, there are N/2 states for k n > and N/2 
states k n < 0. When N is odd, there are (N — l)/2 states 
for k n > and (N — l)/2 states for k n < 0. Additionally, 
there seem to appear two modes with ak n = ±7r at n = 
(N — l)/2. However, they are identified with one another, 
since they are located at the boundary of the Brillouin 
zone. There are N states in both of the cases, as agrees 
with the TBM result @. 

Trigonal symmetry group: The symmetry group of the 
trigonal nanodisk is C^ v , which is generated by the 27r/3 
rotation C3 and the mirror reflection a v . It has the ir- 
reducible representation {Ax, A2, E}. The Ai represen- 
tation is invariant under the rotation C3 and the mirror 
reflection cr v . The A2 representation is invariant under c.3 
and antisymmetric under o~ v . The E representation ac- 
quires ±27r/3 phase shift under the 2-7r/3 rotation. The 
A\ and A2 are 1-dimensional representations (singlets) 
and the E is a 2-dimensitional representation (doublet). 
These properties are summarized in the following char- 
acter table. 





c 


2c 3 


3a v 


Ai 


1 


1 


1 


A 2 


1 


1 


-1 


E 


2 


-1 






(10) 



The mirror symmetry is equivalent to the exchange 
of the K and K' points. With respect to the rotation 



there are three elements c°, C3, c 2 , which correspond to 
1, e 2 ™/ 3 , e 4 ™/ 3 . Accordingly, the phase shift of one edge 
is 0, 27r/3, 47r/3. From this requirement we deduce that 
the state, indexed by the edge momentum k n as in ([9]), is 
grouped according to the representation of the symmetry 
group G% v as follows, 

Ax (singlet) : \k° n ) + \ - k° n ), 1 = 6n + 3^ 
A 2 (singlet): \k° n ) - | - k° n ), / " 3Na ' (n) 

6n i 1 

£ (doublet): |fc±), \-k£), fc rt = ^^^ 
where k" is subject to the condition ([3]). It follows that 
[(N + l)/3\<n<[N/2\, (12) 



where [a\ denotes the maximum integer equal to or 
smaller than a. Some examples read 



N 


3 


4 


5 


6 


7 


Ax 




fc? 




l0 
K 2 


uo 
n 2 


A 2 




fc° 


k\ 


k^ 


h>2 , fi-3 


E 


±k+ 


±fc 2 " 


1 5 1 ^2 


±fe 2 i 


±fc 3 , ifcg" 



(13) 

The numbers of doublets (S-mode) and singlets (Ai- 
mode) are given by [±(N + 1)J and N — 2 [_^(N + 1)J , 
respectively. 

Berry connection: To see the meaning of the wave 
number fc" more in detail, we have calculated the Berry 
connection for various states, which we show for the case 
of N = 7 in FiglHl Each mode is found to be chiral. 
We observe clearly a texture of vortices: The number of 
vortices is 6,7,7,9 for Ik®), \k$), \k$), \k°), respectively. 
The vortex at the center of mass has the winding number 
2 in \k^}. In general, the total winding number N VOI t ex 
is calculated by 

at ~ l I j <t>*(x,y)di(t)(x,y) 

iVvortex = —fdXi — ^ = N + 771- 1, (14) 

2-kJ \<p{x,y)F 

with m = 0, 1,2, ••• ,[(N-l)/2\ in the size-TV nan- 
odisk, where the integration is made along the closed 
edge of a nanodisk. There are n vortices along the y- 
axis in the state \k"). The state being the E- 

mode, has a vortex at the center of mass, where the 
winding number is 2 in the state On the other 

hand, the state \k„) does not have a vortex at the center 
of mass, and the combinations |fc°) ± | — k%) belong to 
the Ai and A 2 representations, respectively. This state- 
ment is demonstrated by investigating the zero points of 
the envelop function ([7]), where vortices appear. For in- 
stance, it is expanded around the center of mass z = zq 
as ip (z) — J2 n =o C n {z — zo) n , where the coefficients Cq 
and C\ are found to vanish at k — k£ and k = k+ , re- 
spectively, with C3 7^ 0. Hence the winding number is 2 
for k = fc+. 

It should be emphasized that there exists a good agree- 
ment with respect to the edge states between the Dirac 
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FIG. 4: (Color on line) The energy spectrum with Coulomb 
interactions in nanodisks with N = 5, 6, derived based on the 
tight-bindind model. The vertical axis stands for the energy 
in unit of leV. The iV-folded degenerate states in the nonin- 
teraction regime split into different levels according to the rep- 
resentation of the trigonal symmetry Ca v , as indicated. For 
instance, the positive (negative) energy levels are for down- 
spin (up-spin) states. The ground state is a ferromagnet. 



description and the exact diagonalization results of the 
TBM even for a small system. Indeed, we can easily 
compute the phase at each lattice point by exact numer- 
ical methods. Then, comparing it with the result due to 
the Dirac description, it is easy to see that a good agree- 
ment holds between them. This shows that the real space 
Barry connection computed by exact numerical methods 
has the same physical reality with the one obtained by 
using the Dirac Hamiltonian. 

Zero-energy splitting due to Coulomb interactions: We 
have constructed explicitly the wave function for each 
state in the size- A trigonal nanodisk. All these states are 
edge modes belonging to the zero-energy sector. When 
Coulomb interactions are taken into account, the degen- 
eracy among the zero-energy states is resolved[6]. The 
Coulomb Hamiltonian has the trigonal symmetry Ca v , 
and the energy spectrum splits into different levels ac- 
cording to its representation, as illustrated in Fig0J 

There exists additionally the spin degeneracy in the 
noninteracting Hamiltonian: The total degeneracy is 2N. 
The spin degeneracy is broken spontaneously due to 
the exchange interaction when Coulomb interactions are 
introduced j6j. The splitting is symmetric with respect to 
the zero-energy level. At half-filling, electrons with the 
identical spin fill all energy levels under the Fermi energy. 
Then, the spin of the ground state is N/2, and it is a fer- 
romagnet. We show the energy spectrum for N = 5, 6 in 
FigJU By tuning the chemical potential any of them is 
made the ground state. 

Magnetic vortices: The zero-energy degeneracy is re- 
solved by the Coulomb interaction, and the dispersion 
relation becomes nontrivial. The time-dependent solu- 
tion is well known, 



^f{t,x) =e 



— i(E"t — k^x) — k^y 



D — i(E^t-\-k™x) — k™y 



(15a) 
(15b) 



where E" is the energy of the state \k") [Fig|3]. Here 
we have suppressed the spinor part. On one hand, the 
Aj-modes |fc°) ± | — fc") represent standing waves. On the 
other hand, the i?-modes \k@) and | — k^) represent the 



right-propagating mode and the left-propagating mode, 
respectively, for /3 = ±. 

Charged particles propagating along a closed edge gen- 
erates magnetic field. The electromagnetic interaction is 
described in terms of the electromagnetic potential A, 
which is introduced to the system by way of the Peierls 
substitution dj —> dj + ieAj/h. From the Weyl equation 
([3]) we derive 



eAi(x) = hAi(x)/\4%(x)\* 



(16) 



with Ai(x,y) = —i(f>"&* (x)di<fr£(x) in the lowest order of 
approximation, where </>a( x ) * s assumed to be not modi- 
fied from ([7]). The potential Ai(x) exhibits the same tex- 
ture of vortices as in FigGU The magnetic field is given 
by 

B(x) = V x A(as) = — V v n 6 (z - z n ) , (17) 
e '-^ 

n 

where v n stands for winding number of the vortex at z = 
z n . Hence a texture of vortices in the Berry connection 
leads to a texture of magnetic vortices. A comment is 
in order. This 6- function type magnetic field would be 
smoothed out in a rigorous analysis of the coupled system 
of the Maxwell equation and the Weyl equation. 

It is intriguing that, by tuning the chemical potential, a 
vortex with the winding number 2 emerges in the ground 
state As is well known, a single flux quantum has 

experimentally been observed in superconductor by us- 
ing an electron- holographic interferometry[16]. Then, in 
principle it is possible to observe a vortex texture in nan- 
odisk as well. Furthermore, by attaching a superconduc- 
tor film one may observe a disintegration of a vortex into 
two when the flux enters into the superconductor from 
the nanodisk. This would verify the winding number 2 
of a vortex. 

Conclusions: In this paper we have classified the zero- 
energy sector of the trigonal zigzag nanodisk into a fine 
structure according to the trigonal symmetry group C$ v . 
We have explicitly constructed wave functions based on 
the Dirac theory and specified them by the quantized 
edge momentum. A texture of magnetic vortices is found 
to emerge, which has an unusual winding number. As far 
as we are aware of, the vortex with the winding number 
2 has never been found in all branches of physics. This 
is because two vortices with the winding number 1 have 
lower energy than one vortex with the winding number 
2 in general. In the present case the disintegration of a 
vortex into two is prohibited by the trigonal symmetry. 
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